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Abstract 

For the past three years, the MIT Haystack Observatory and the broadband team have been de- 
veloping a proof-of-concept broadband geodetic VLBI microwave (2-12 GHz) receiver. Also on-going 
at Haystack is the development of post-correlation processing needed to extract the geodetic observ- 
ables. Using this processing, the first fully-pliase-calibrated geodetic fringes have been produced from 
observations conducted with the proof-of-concept system. The results we present show that the phase- 
calibrated phase residuals from four 512 MHz bands spanning 2 GHz have an RMS phase variation of 
8° which corresponds to a delay uncertainty of 12 ps. 


1. Raw Fringe Phasor Model 

This development begins by assuming that one has at their disposal the normalized correlation 
coefficients for each frequency channel (whether represented in the cross-power frequency or lag 
domain) for each single correlator accumulation period (AP). The raw fringe phasor <& r (/) is then 
defined as the phase of the cross-power spectrum produced from the correlator output data and is 
assumed to be the product of the fringe phasor d> 3 (/) due to the residual geodetic delay and the 
hardware (phase cal) phasor & pc (f), where / is the sky frequency. Furthermore, since amplitudes 
are currently not incorporated in the phase cal processing, the magnitude of all phasors considered 
in this note are arbitrarily set to unity. 

The hardware-related phase is given by the phase calibration phasor $ pc (/) which represents 
the arithmetic difference between the phase cal phases obtained at each station. <h pc (/) is also 
sampled every / = 5n MHz, where n is the rail harmonic number. In the context of this note, 
% c (f) will be modeled as a piecewise-linear phasor function composed of as many pieces as there 
are frequency channels as indicated in Figure 1. It is important to note that d> pc (/) is not simply 
the difference of the phase cal phases obtained directly from the correlator, since the peal phases 
are sampled every 5 MHz and the fringe phases in each frequency channel, generally, are sampled 
at a different rate. Additional processing (interpolation) is needed to obtain <I? pc (/) from the 
correlator peal phases. This is discussed in the section on phase calibration fitting. 

In regards to the phase due to the residual geodetic delay, systematic errors such as those 
introduced by clock drift and ionosphere delay are neglected here to place emphasis on the phase 
calibration processing needed to connect the raw fringe phases across several receiver bands. The 
residual geodetic fringe phasor (i.e., the error-free fringe phase), represented by the function <F p (/), 
is assumed to be linear in frequency since the aforementioned errors are ignored; estimation of d> 9 (/) 
from knowledge of d> r (/) and <3? pc (/) is the goal of the correlator post-processing described in this 
paper. $ 9 (/) can be formally written as: 

%(f) = e 2 " /Ts (1) 
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where r g is the residual time delay (between the two stations) to the geodetic model. Under this 
development the raw fringe phase is assumed to be the product of the geodetic and phase cal 
phasors: 

Mf) = <M/)<M/) ( 2 ) 

Figure 1 displays a graphical example of the phases d v(/), < h pc (/), and 3> s (/) for the current 
receiver architecture (eight 32 MHz channels per band and four 512 MHz bands) and band-to-band 
spacings Afi through A/ 3 . 



Figure 1. Graphical representation of geodetic fringe, raw fringe, and phase calibration phases. 


2. Phase Calibration Phase Fitting - Derivation of tpc(F) 


As mentioned previously, the phase cal tones are injected in the front-end of the receiver at 
frequencies which are multiples of 5 MHz whereas the frequencies at which the raw fringe phases 
are sampled, in general, are different. As a result, the phase cal data provided by the correlator 
must be fit to a model (linear in this case) in order that the phase calibration phases can be 
evaluated at the same sky frequencies as the raw fringe phases. The following model is used to 
describe the phase cal phasor provided by the correlator: 


<| ) m (f) _ e 2nnf r T™+'t>™ 

P f=nfr 


(3) 


<&£(/) describes the phase cal model for the mth frequency channel as a function of tone index 
n; in this model <3?™(/) represents the arithmetic difference of the phase cal phases (as provided 
by the correlator) between the two stations for the given channel. In equation (3), f r is the rail 
frequency spacing (5 MHz) and r™ and are the phase cal delay and phase, respectively, for 
the mth channel; / and therefore n are bounded by the frequencies of the channel from which the 
phase calibration tones were extracted. Given the samples of <h™(/), the task is to estimate r™ 
and which, by inspection of equation (3), can be accomplished by applying a discrete Fourier 
transform (DFT) to <I?™(/). The DFT makes no assumptions about the frequency locations of 
the tones in any given channel, so there is no restriction on the tone selection. After r™ and 
4> 7 p C have been estimated, the interpolated peal phasor d>™(/) can be expressed using the same 
linear model in equation (3) by simply substituting for the tone frequency nf r the sky frequency 
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/ corresponding to that at which the raw fringe phasor is evaluated: 

$™(/) = (4) 

The complete phase calibration phasor (assuming there is no overlap of the frequency channels) is 
just the sum of the interpolated peal phasors for each frequency channel: 

<M/> = E 3, ?(/) d) 


3. Phase Calibration Correction of Raw Fringe Phases 

The raw fringe phasor 4> r (/) and phase calibration phasor <f>p C (/) are produced by the correlator 
after each accumulation period (AP); therefore, <f>™(/) can be corrected for each AP independently 
(this is the so-called AP-by-AP peal mode). The benefit of the AP-by-AP mode is that hardware 
drift over extended scan periods is removed. Equation (2) indicates that the residual geodetic 
phasor can be obtained by simply multiplying the raw fringe phasor by the complex conjugate of 
the phase calibration phasor: 

<W) = MfWpM) (6) 

For the polyphase filter implementation of the frequency channelization currently implemented in 
the digital backend, the sky frequency / can be decomposed into a summation of the LO frequencies 
of the channels and the cross-power frequency relative to the channel LO frequency. Formally put: 

/ = ffi) - kAf cpf (7) 

where f™Q is the LO frequency of channel m, f cp f is the cross-power frequency sample interval, 
k is the cross-power frequency index (i.e., 0,1,2,...A SC -1), and N sc is the number of cross-power 
spectral samples in a single 32 MHz frequency channel. Equation (6) can now be rewritten using 
equations (4-5,7) to provide a more detailed formulation of the method used to correct the raw 
fringe phases and obtain the residual geodetic phases: 

$ g (f) = «M/) e j2TrkAf cpf T™ e -j4 (g) 

m 

4. Broadband Phase Calibration Processing of Experiment 3296 

Since the goal of the broadband phase calibration processing is to remove the hardware-related 
delay and phase from the raw fringe phasor, broadband fringe tests have been conducted on the 
Westford-GGAO baseline to demonstrate the performance of this processing. The broadband 
system possesses the flexibility to change the inter-band spacings Af\ through Afz (depicted in 
Figure 1) with 400 kHz precision. In fact, the inter-band spacings can be set identically to zero 
which provides a quite useful diagnostic for validating the broadband phase calibration process- 
ing as outlined in this note. In this so-called “overlapping bands” configuration, all four bands 
observe an identical set of frequency channels, and the residual geodetic fringe phasors for each 
band (fringe fit independently) are expected to be identical if the receiver is truly front-end noise 
limited. In practice, the residual geodetic fringe phasors must be coherently averaged over all 
APs in the integration period in order to detect the interference fringes. Table 1 displays the raw 
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Table 1. Raw singleband/multiband delays, fringe rate, and phase results from fringe fitting. 



A 

B 

C 

D 

Singleband Delay (ns) 

-2.334 

1.374 

-10.086 

-4.063 

Multiband Delay (ns) 

-1.970 

1.768 

6.011 

-3.635 

Fringe Rate (mHz) 

1.522 

2.216 

2.016 

2.043 

Fringe Phase (°) 

-157.7 

-91.5 

55.6 

-72.2 


Table 2. Phase-calibrated singleband/multiband delays, fringe rate, and phase results from fringe fitting. 



A 

B 

C 

D 

Singleband Delay (ns) 

-0.503 

0.180 

-0.223 

-0.241 

Multiband Delay (ns) 

0.076 

0.069 

0.083 

0.092 

Fringe Rate (mHz) 

1.578 

1.576 

1.586 

1.588 

Fringe Phase (°) 

-144.5 

-142.5 

-144.6 

-143.6 


singleband/multiband delays, fringe rate, and phase for HH polarization obtained from an obser- 
vation conducted on Day258 at 1403UT year 2009 (Haystack Exp: 3296) coherently averaged over 
the entire 10 minute scan. In this experiment, all four bands observed 6408.4-6888.4 MHz. The 
raw fringe phasors processed without the broadband phase calibration method are given in Table 
1 and the calibrated counterparts of those parameters are given in Table 2. The uncertainties 
associated with the singleband/multiband delay, fringe rate, and phase are 0.2 ns, 0.013 ns, 0.011 
mHz, and 1.4°, respectively. Comparison of Tables 1 and 2 demonstrates that the phase calibration 
processing is performing nearly as expected since the calibrated results for all four bands are in 
such good agreement; the minor discrepancies are under investigation. 

During this same experiment, data were also collected such that the four bands spanned a 
total contiguous bandwidth of 2 GHz (6.4 — 8.4 GHz). Differences in cable lengths and local 
oscillator phases in the hardware comprising each frequency band introduce a hardware-related 
component to the raw residual fringe phase shown in Figure 2a. The discontinuities in Figure 2a 
appear at the spectral boundaries of adjacent frequency bands as a result of these hardware-related 
components. The raw data were phase calibrated using the technique described above, and the 
resultant calibrated fringe phases are shown in Figure 2b. As is obvious in comparing Figure 2a to 
2b, the discontinuities are suppressed in the calibrated geodetic fringe phases. The resultant mis 
phase noise is 8.5°, which translates to a delay uncertainty of 12 ps across the 2 GHz of contiguous 
bandwidth. This delay uncertainty will be reduced by observing larger spanned bandwidths and 
this is one of the explicit goals of VLBI2010. 

5. Conclusions and Future Plans 

In order to remove the hardware-related phase and delay components from the raw fringe 
phasors, a new phase calibration correction algorithm has been developed. This processing was 
shown to reconcile the band-to-band delay and fringe rate differences when all four bands observed 
the identical set of sky frequencies, thus demonstrating the algorithm’s ability to remove the 
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Figure 2. Raw and Calibrated Residual Fringe Phase vs. Sky Frequency 


hardware-related components of the aformentioned delay/rate parameters. Though the observable 
parameters are fairly consistent, Table 2 demonstrates that there are post-correction differences 
between the bands. Observation of such discrepancies suggests that additional independent noise is 
added to the signal recorded in each frequency band after the 2-12 GHz microwave signal received 
by the front-end is equally split into each of the four frequency bands (signal chain [1] pp. 35). 
A signal chain cascade analysis indicates that the noise contributed by the analog components in 
each frequency band is negligible. The noise contributed by the digital backend (DBE) (i.e. , ADC, 
fixed point math quantization errors), however, is not currently a quantity that can be measured 
directly. As such, there may be significant noise contribution by the DBE of which we are unaware. 
This situation will be rectified with the introduction of the next-generation DBE, at which point a 
full signal chain (analog/digital) noise budget can be developed to assess the expected uncertainties 
in the band-to-band post-correction residuals. 

The phase calibration correction algorithm was also shown to align the raw fringe phases across 
a total bandwidth of 2 GHz, which resulted in a delay uncertainty of 12 ps. In order to reduce 
the delay precision of the geodetic observables, future observations will be conducted which span 
bandwidths greater than 2 GHz. In observing these wider bandwidths the ionosphere is expected 
to introduce a dispersive (non-linear) behavior in the fringe phase function (e.g., Figure 2b) which 
will require modification of the fringe fitting routine from a linear model to one that is dispersive. 
Finally, the results reported in this paper were for the linear HH polarization only. Developments 
of the fringe fitting algorithm will also need to incorporate all four polarization components (HH, 
HV, VH, and VV) into the geodetic fringe fitting process in order to further enhance the SNR and 
realize the goals of VLBI2010 [1]. 
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